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SUMMARY 


Using the nonlinear load elements in thermal analysis to simulate an 
undocumented nonlinear thermal boundary condition is presented. The treatment 
of the nonlinearity arising from the temperature-dependent convective film 
coefficients is shown in detail. As an illustration, emphasis is placed on 
the modeling techniques and their interrelationships with the solution accu- 
racy as affected by a specific integration algorithm of the transient thermal 
analysis used in the NASTRAN Thermal Analyzer. Briefly shown is the under- 
lying theory on which the maneuvering of terms pertinent to the modeling 
depends. This demonstration provides some insight into the intricacies of the 
method that would be general to all applications. A recommendation is also 
made to modify a nonlinear load element that will enhance the solution capa- 
bility and broaden the scope of application. ' 


INTRODUCTION 


The known solution capabilities of the NASTRAN Thermal Analyzer (NTA) 
have been well documented (ref. 1 and 2), and a large part of them have been 
demonstrated in detail (ref. 3). Despite their frequent encounterment in 
engineering applications, certain nonlinear solution capabilities have not 
been provided by the NTA in the transient-state thermal analysis (APP HEAT, 

SOL 9). Those excluded are, for instance, the temperature-dependent thermo- 
physical properties, temperature-dependent convective film coefficients, etc. 
Their absence from the original NTA's capability list was solely attributed to 
the fact that the required computer processing time (CPU and 1-0 times) would 
be excessive. Computations involving the reassemblages of the thermal conduct- 
ance matrix, together with the decompositions of the relevant matrices that 
contain the nonlinear thermo physical properties, would have to be repeated for 
all steps of the time increment in integration. The program shortcomings, 
however, may be compensated by an appropriate application of the available 
nonlinear load elements (NLLEs) in NASTRAN, although some elements were origi- 
nally developed for structural analyses. The intrinsic modular structure of 
NASTRAN permits them to be accessible to the NTA. Those NLLEs, obviously, can 
function as nonlinear thermal loads, and they can be employed to simulate other 
complex physical phenomena as well. Otherwise, new modules would have to be 
developed and implemented, even though an excessive computer processing time 
was tolerable. Such an alternative would be an expensive and time-consuming 
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proposition. The user would be deprived of immediate solution meanwhile. 

Six NLLEs are currently available in NASTRAN, namely N0LIN1*, N0LIN2, 
N0LIN3, N0LIN4, N0LIN5 and NFTUBE. The last two were developed specifically 
for thermal analyses: N0LIN5 is capable of simulating temperature sensitive 

thermal radiation surface properties, such as a louver (ref. 4), and NFTUBE is 
capable of dealing with thermal energy transfer between the wall and the fluid 
flowing inside a tube (ref. 5). The remainders are four nonlinear elements of 
general type capable of accommodating four different functions with a temper- 
ature or heat flux being the argument. N0LIN1 allows to specify a table through 
TABLED! (i = 1,2, 3, 4), and it is the most versatile NLLE in NASTRAN. 

The main purpose of this paper is to illustrate how to acquire a non- 
routine, but essential solution capability by an application of the NLLEs 
through modeling. For a systematic demonstration, N0LIN1 was employed to 
simulate the temperature-dependent convective film coefficients, h(T). The 
importance of this problem is evident in view of a variety of applications 
ranging from a domestic heating system, a hot gas flowing through a nozzle, 
a fluid flow in a nuclear reactor, and a flow in a wind tunnel during the 
start-up to a thermal ablation including a receding surface. They are all'" 
characterized by such a nonlinear boundary condition. More specifically, 
treatments of the temperature-dependent film coefficients based on one and 
two reference temperatures are given in detail. The dilemma arising from the 
presence of the MPC to establish an average temperature of two reference 
temperatures and the reference with that average temperature by the NLLEs was 
resolved by a trick modeling. A brief review of the theoretical aspect serves 
not only to provide a basis for the maneuvering in modeling and to provide 
insight into the intricacies of the method, but to amplify basis of modeling 
techniques that are equally valid to other applications. A recommendation for 
the improvement of the N0LIN2 is outlined. It would enhance solution capa- 
bilities and broaden the scope of engineering application. Guidelines for the 
user to apply the NLLEs in thermal analysis are provided. 


THE DEMONSTRATIVE PROBLEM 


For the convenience of easy reference and discussion, and later on, to 
demonstrate the modeling techniques, a relatively simple problem without losing 
any generality has been designed, and it is described below: 

The transient temperature responses in a composite slab (pipe) of an infinite 
extent are to be determined. Initially, the slab has a uniform temperature of 
293 K (20 C) throughout its body. The boundary surface of the metallic layer 
is in contact with a stream of hot gas maintained at 1422°K (1149°C) from time 
t = 0. The dimensions and thermophysical properties of the slab are tabulated 
in Table 1. The system depiction, together with its finite element model, 
is shown in Figure 1. 


*The names of actual NTA cards are capitalized. 
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THE THEORETICAL CONSIDERATION 


The success of attempting a solution utilizing any NLLEs to simulate a problem 
with nonlinear thermophysical properties relies upon the following factors: 

(1) The nonlinear terms, arising from the temperature-dependent thermophysical 
properties, which are embraced in the square matrices of the thermal conduct- 
ance and/or the thermal capacitance, must be extractable and separable from 
their normal positions. The resulting products of the temperature-dependent 
icoef f icients (kA/£ an d hA or C) and their respectively associated variables 

i(T or t) must be transposed and merged into the nonlinear thermal load . 

vector . 

(2) An NLLE of an appropriate form suited to accommodate the transposed 
quantities must be selected. Any restraints associated with the use of the 
NLLEs must be strictly observed. 

(3) A full acquaintance with the characteristics of each NLLE and the 
integration algorithm for the transient heat transfer analysis (APP HEAT, SOL 9) 
is essential in controlling and accomplishing desired solution accuracy, 
stability and efficiency at wish. 

To serve as the theoretical basis for such a maneuvering, a brief review 
of pertinent mathematical operations is in order. 

A transient thermal analysis in the NTA is based on the general heat equation 
in the matrix form as follows (ref. 1): 

[ C ] { T } +[K]{T} = {Q f } + {Q n } (1) 

where 

T = the vector of temperatures at grid points 
T = the rate change of temperature vector 
K = the thermal conductance matrix of constant elements 
C = the thermal capacitance matrix of constant elements 
= the vector of applied linear thermal loads 
Q n = the vector of applied nonlinear thermal loads 

The convective coupling terms are normally incorporated with the thermal 
conductive couplings in the thermal conductance matrix [K] regardless of their 
dependency. When the convective film coefficient is temperature-dependent, the 
convective coupling terms embracing h(T) have to be separated from their orig- 
inal [ K„ J , and thus[K„] |Tj. The extracted terms being a new vector |h(T)| 
are transposed to the other side of the equal sign in equation (1) , and they 
become or merge with the vector of nonlinear thermal loads, {Q n }* The 
rearranged new equation conforms exactly to equation (1), meaning that it is 
then within the solution capability of the NTA. , Mathematically, it can be 
expressed by 
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( 2 ) 


[K 0 (k,h(T))][T} = [ K-l (k) ] {t } +{h (T)| 

The integration algorithm for equation (1) is based on the modified 
Newmark Beta Method (ref. 6). It is a forward and backward differencing 
method with a parameter 0 that enables the user to select a value in the 
range of O<0<1. When equation (1) is transformed into such a difference 
expression, it reads (ref. 1): 

(K]{dT„+l + <l-f»T ll }+ JJC]{t „ +1 - T n } 

- {uq ' +1 + a-tn {qA>} + a+» K}-f {Cl} «> 


the subscript n refers to the nth time step. The accuracy and stability of 
solution are closely related to 0 . Equation (3) can be rearranged to give 

[l(C] + fi[K]]{T n+1 }=[i[C] - (l-0)[K]]{T. n | 

At L At 

+ {eQn+l + (l-0)Qn} + a+0){Q£} -0{Qn-l} < 4 ) 

The mat ri x [1 / At C + 0 K] is d e c om posed i nto its t riangular^ factors^ and 
solved at each time step via the forward and backward substitution 
process. If [K] contains terms of k(T) or h(T), a search from the tables 
and replacement of new values based on {T n } have to be carried out for all 
terms residing inside those two square matrices that are associated with 
|T n +jJ and {T n |. Consequently, the matrix decomposition has to be performed 
for each time step. This requires a very excessive computer processing time. 
However, [K] needs be evaluated only once if its constituent terms are all 
constant. This situation is equally valid in [C] , although C(T) is not 
considered in this study. 

Relating the terms of nonlinear thermal load to the specific demonstrative 
problem, it is obtained from equations (1) and (2) that 


(5) 


( 6 ) 


{q 11 } = - [h(T>] 


J \ 


where 


0 
0 

' H 200' 


{ 


H x = h(T r )AT L 

h 200 = “ h( T r) AT 200 
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and h(T r ) denotes the dependency of the convective film coefficient in a refer- 
ence temperature T r which can be the wall temperature or an average temperature 
of the free stream and the wall temperature as will be discussed and treated in 
following sections. In the NTA, the N0LIN1 card defines an NLLE of the form 

q"(t) = s^Cr-jO:)) (7) 

where Q^(t) is the load being applied to the GRID i at time t, S-^ is a scale 
factor. F(Tj (t) denotes a tabulated function in Tj to accommodate H^(Tj) 
and ^(Tj) as defined in equation (6). The equivalent nonlinear thermal load 
functions are entered through a TABLEDi (i = 1,2, 3, 4) card. T (t) is any 
permissible time-varying temperature on which the functions depend. Other 
NLLEs, N0LIN2, N0LIN3 and N0LIN4, can be selected as long as their specific 
forms are suitable for respective applications. 


THE CONVECTIVE FILM COEFFICIENT h(T r ) 


The convective film coefficient h governs the magnitude of heat transfer 
between the hot gas and the boundary surface of structure, h, in general, can 
be expressed as h(T r ) where the reference temperature T r may refer to only one 
temperature, such as the wall temperature, or more than one temperature through 
an intermediary, such as the average value of the free stream temperature and 
the surface temperature of the wall. A typical example is a fully developed 
pipe flow. The expression for h is (ref. 7) 

h = 0.023 hz- (Re)°' 8 (Pr)°' 4 (8) 

Dh 


For large temperature differences, a mean temperature is used in evaluating 
all physical properties except the flow velocity. This procedure would involve 
very tedious interrations. To simplify computations, it is desirable to 
separate the temperature-dependent from the independent factors. An alterna- 
tive form is, therefore, obtained by an application of the equation of state 
and the continuity equation, which leads to 


h 


= 0.023 C r 


0. 2 - 0.6 0.8 - 0 . 2 . 0.8 

./!» Pr« G» Dh 0F> 

J- oo 


(9) 


As the free stream temperature T B is known, all quantities in the preceding 
equation except T w are constants, h is, therefore, a function of T w . For the 
problem of concern, it is h = h(T^). The values of this function obtained by 
using equation (9) are listed in Table 2. Similar expression, but referencing 
with an average value of the free stream temperature and the surface tempera- 
ture is also widely used in engineering practices. For this demonstrative 
problem, it is h = h(T r ) and T r = %(Tj + T 200)- The values of h vs T r are 
given in Table 3. 
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THE MODELING TECHNIQUES AND DISCUSSION 


A one-dimensional finite element model is sufficient to represent the 
demonstrative problem, and this model is shown in Figure 1, A listing of the 
input data deck is reproduced in Figure 2. With comments appropriately inserted, 
functions regarding various segments of the deck are self-explanatory. The 
space bounded by two rows of asterisks is reserved for adding a packet of cards 
which simulate the effect of the convective film coefficients. Three cases of 
h are considered: (1) h of a constant value, (2) h = h(Ti) , and (3) h =h(T r ) 

with T r = h(Ti + T200)- 

h of a constant Value 


To assure correctness of modeling and to evaluate the obtainable accuracy 
of solution by the nonlinear load approach, the case of a constant film co- 
efficient with h = 0.467 W/cm^-K was tested first. Because, this case is well 
within the documented capability of the NTA. It can be accomplished using a 
train of cards CHBDY, PHBDY and MAT4 by the conventional modeling method. The 
temperature results so obtained serve as a base to compare with that obtained 
by the nonlinear load approach. The two packets of input data cards, whose 
images are shown in Figures 3(a) and 3(b), were used in respective computer 
runs. The equivalent nonlinear thermal loads for the pair of NLLEs through 
N0LIN1 405 with TABLED1 4051 and TABLED1 4052 are given in Table 4. The effect 
of selecting the time-step size At for integration on solution accuracy was 
also studied. 

Temperature results summarized in Table 5 with the headings of Cases 1 and 
3 are solutions yielded by these two approaches. In both cases, the same time- 
step size (the At-set (a) in Figure 4) was used. The largest difference between 
the two temperature profiles occurs at t = 0.25 sec., suggesting that a refine- 
ment of the time-step size is warranted. To maximize the efficiency of the 
computer processing time, refinements of At were made only in the beginning 
segments of time. Three sets of time-step size, which are presented in the 
NTA TSTEP card format and identified as the At-sets (a), (b) and (c) , are shown 
in Figure 4. A total of five cases were computed, two (Cases 1 and 2) by the 
conventional modeling method and the other three (Cases 3, 10 and 11) by the 
nonlinear load approach. Temperature results are tabulated in Table 5 for 
comparison. It is evident that both modeling methods can yield comparable 
temperature solutions with a desired accuracy, provided an appropriate At-set 
is used (e.g. , to compare Case 2 with Case 11). 

This test gave confidence in the nonlinear load approach and verified the 
adequacy of the two sets of time-step size, namely the At-sets (b) and (c) . 

They were, therefore, selected in all following computations. Inevitably, 
efficiencies of those runs using the At-sets (b) and (C) were penalized by a 
certain amount. They are reflected by the CPU and 1-0 times, which are also 
summarized in Table 5. 
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h = h(TT) 


The modeling procedure similar to that used in the preceding case of a 
constant h is equally applicable to this case of a convective film coefficient 
h, which is dependent on the surface temperature of the wall. The preparation 
of the equivalent nonlinear thermal loads, however, is much more involved. 
Computations of those equivalent loads were based on the expressions given in 
equation (6). The computational steps and results are detailed in Table 3. 

The relevant packet of input data cards is shown in Figure 3(c). 

Two computer runs were executed with the specified At-sets (b) and (c) . 
Table 6 displays their temperature results under the headings of Cases 12 and 13. 
Results are in excellent agreement, meaning that either At-set is satisfactory 
to render accurate solutions. The At-set (b) was then selected in studying the 
effect of the parameter P in the integration algorithm on temperature solutions. 

The integration algorithm available in the NTA for its transient thermal 
analysis is a modified Newmark Beta Method. It uses a fixed time-step size, 
but permits different time-step sizes to be changed at discrete times upon the 
user's specification. To avoid excessive matrix decompositions, the number of 
changes should be kept low. This algorithm is known to be an efficient one for 
large systems of linear equations, but it is uncertain for problems with strong 
nonlinearities. As accuracy and stability of a solution are closely related to 
the parameter (3 , a numerical investigation with different values of P was 
conducted. The values of p , including 0, 0.25, 1 and 0.55 (the default value), 
were executed. Except the case of using p = 0.5 (the central differences), 
whose temperature results are virtually identical to that of Case 13 (using p 
= 0.55) and, therefore, eliminated from the tabulation in Table 6. Case 16 
(using p = 1.0, the backward differences), yields comparable temperature results 
as that of Case 13. However, the solution instabilities caused both cases with 
P =0.25 and 0 (Euler integration) to terminate in their mid-courses. They 
are respectively labeled as Cases 17 and 18 and also included in Table 6. 

h = h(T r ) with T r = h(^l + T 200 ) 

When the average temperatures T r of the free stream temperature T 20 O anc * 
the surface temperatures of wall Ti are referenced by the temperature-dependent 
convective film coefficients h = h(Tr) , the MPC (Multipoint constraint) card is 
suited to provide such a relation. When T 100 = Tr i s assigned, the equivalent 
nonlinear thermal loads are functions of Ti00> meaning that Tioo(t) will be 
referenced by the nonlinear thermal load tables, if the modeling technique 
similar to that used in the two preceding cases is followed. Since Ti00 is 
constrained by an MPC, this modeling has violated one cardinal rule that must 
be strictly observed by a user in applying the nonlinear loads. The grid points 
to which the nonlinear loads are applied and the variables (either temperature 
or heat flux) on which they depend must not be variables eliminated by con- 
straints, i.e., they must be in the solution set. To circumvent this dilemma, 
a fictitious point GRID 150 was created for T 150 , and it was connected to T100 
by a super conductor, so that Tioo(t) = Ti5o(t). T 150 instead of T 100 was then 
referenced by. h(T r ), i.e., h(T r ) = h(Tioo) = MT 150 ), which, in turn, produced 
the equivalent nonlinear thermal loads. Figure 3(d) shows all relevant cards 
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and their relationships in detail, and the entered quantities in the nonlinear 
load tables are easily identifiable in comparison with those values given in 
Table 4. 

The two t-sets (b) and (c) were again selected to execute the computer 
runs. Results with the labeling of Cases 14 and 15 are shown in Table 7. 
Temperatures are still overall in good agreement with each other, but they no 
longer show the same quality as evidenced in previous cases. A scrutiny of the 
results reveals that relatively large temperature differences exist between 
these two temperature profiles in the range of time from 1.0 to 3.0 sec. This 
is a typical example showing the intricacies of this method , the interplays of 
the selected time-step sizes with the integration algorithm, and the effect on 
temperature results due to the presence of the MPC. 

Equation (4) shows that an inherent time lag error exists in applying the 
nonlinear thermal loads. They are computed for T of the (n+l)th time step 
based on temperature values at the nth and (n-l)th time steps. The larger the 
time-step size, the greater the deviation of the nonlinear thermal load from 
the desired value. When MPC is employed to relate a floating temperature (e.g., ; 
T]^) to a reference temperature (e.g., T 150 ) on which values of the nonlinear 
thermal loads depend indirectly, it introduces an additional step of time lag, 
and thus, an aggregated error for the nonlinear thermal load value. The errors 
are now induced by the time lag steps spanning from the (n+l)th to (n-2)th 
steps. Consequently, the selection of a finer time step size is appropriate 
in applying the nonlinear thermal loads, especially when an MPC is present 
to average two temperatures. 


OTHER APPLICATIONS AND RECOMMENDATIONS 


The modeling techniques of using the NLLEs to simulate the temperature- 
dependent convective film coefficients have been methodically illustrated. 

The demonstrated equivalent nonlinear thermal loads are all functions of 
temperature. If a simulation calls for the nonlinear thermal load to be a 
function of heat flux, a combination of the EPOINT and TF cards, in addition 
to the NOLINls and TABLEDls, have to be used. In all cases, values for the 
equivalent nonlinear thermal loads must be prepared by the user manually. This 
is an error-prone and elaborate process inherent in using the NLLEs. This 
deficiency, however, can be removed with only a relatively simple change made 
to the N0LIN1. It -will further broaden its capability greatly. The recommen- 
dation is to modify the computer code pertinent to S, the fifth field, in the 
N0LIN1 card so that it will accept any time- varying temperature, whether 
specified or computed values at some grid points in the model* The expanded 
capabilities will include, for instance, the same demonstrative problem with 
the hot-gas temperatures being a specified time-dependent temperature function. 
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CONCLUSIONS 


The nonlinear load elements have been illustrated to simulate nonlinear 
boundary conditions, in particular the temperature-dependent convective film 
coefficients in transient thermal analyses. The following remarks also serve 
as guidelines for those who attempt to use the NLLEs: 

(1) Prior to applying the nonlinear thermal loads to an NTA model in a 
transient thermal problem, a trial run on a simplified model is advisable. 

Only the linear thermal loads equal in magnitude to that of the nonlinear load 
set are necessary to apply to this model, and a relatively coarse time-step 
size may be used in the execution. The profile of the temperature results will 
provide a good indication on which the final selection of appropriate time-step 
sizes in runs with the real NLLEs depends. 

(2) Taking into consideration of the time lag errors, very small time-step 
sizes are always desirable. However, it is not necessary to select the same 
fine size uniformly over the entire range of time for integration. To maximize 
the processing efficiency, the selection of nonuniform time-step sizes according 
to the results of the trial run is a sensible and practical approach. Very fine 
time-step sizes should be reserved to concentrate at which the temperature 
profile has the largest slope in its rate change of temperature curve. 

(3) When MPC (also TF, if used) is present in the model for establishing a 
reference temperature that is referenced by the NLLEs, such as N0LIN1, a ficti- 
tious super thermal conductor has to be used to connect to that original grid 
point of referencing. A new grid point at the other end of the super thermal 
conductor is then to be referenced by any NLLEs so as not to violate the cardi- 
nal rule in applying the nonlinear loads. 

(4) For solution stability, which also affects accuracy, the proper selection 
of a value for the parameter (3 of the integration algorithm depends on the 
degree of nonlinearity in presence. However, only the range of 0 . 5 < |3 < 1 
should be considered. 
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TABLE 1 

DIMENSIONS AND THERMOPHYSICAL 
PROPERTIES OF THE COMPOSITE SLAB (PIPE) 


Material 

Stainless Steel (a) 

Insulation- Phenolic (b) 

Thermal Conductivity 
k (W/ cm-k) 

0.16258 

0.00502 

Thermal Capacitance 
PCp (J/cm^-k) 

4.029 

1.725 

Thickness 
Ax (cm) 

0.30 

0.40 

Cross-section 
A (cm^) 

1.0 

1.0 
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TABLE 2 


THE EQUIVALENT NONLINEAR LOADS FOR 
THE CASE OF h = 0.467 W/cm 2 -k 


T200 

N200 = 

— 

hAT200 

_ 

Ti 

-N]_ = hATi 

1422.0 

664. 

074 

293 

136.831 





589 

275.063 





867 

404.889 





1144 

534.248 





1422 

664.074 


TABLE 3 

THE EQUIVALENT NONLINEAR LOADS 
FOR THE CASE OF h = h(Ti) 


Tl 

h(Ti) 

-Nl(Ti) = hATx 

t 200 

N 200( T l) = hAT200 

293 

0.132 

38.676 

1422 

187.704 

589 

0.231 

136.059 



328.482 

867 

0.314 

272.238 



446.508 

1144 

0.393 

449.592 



558.846 

1422 

0.467 

664.074 



664.074 
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TABLE 4 


THE EQUIVALENT NONLINEAR LOADS 
FOR THE CASE OF h = h(T ave ) 


t 200 

Tl 

T ave 

k( T ave-* 

N 200( T ave ) = hAT 200 

- N l< T ave> " hAT l 

14: 

12 

293 

857.5 

0.132 

187.704 

113.190 



589 

1005.5 

0.231 

328.482 

136.059 



867 

1144.5 

0.314 

446.508 

272.238 



1144 

1283.0. 

0.393 

.558.846 

449.592 



1422 

1422.0 

0.467 

664.074 

664.074 













TABLE 5 


A COMPARISON OF TEMPERATURE RESULTS 
FOR THE CASES WITH A CONSTANT h 



Case No. I 

Time 

(see) 

1 

2 

3 

10 

11 

0 

293.00 

293.00 

293.00 

293.00 

293.00 

0.25 

575.82 

580.41 

578.73 

582.08 

581.00 

0.50 

665.64 

668.29 

667.90 

669.82 

668.75 

0.75 

723.79 

725.67 

725.68 

727.06 

726.06 

1.00 

768.38 

769.88 

770.07 

771.19 

770.23 

1.25 

806.14 

807.42 

808.22 

809.19 

807.75 

1.50 

839.92 

841.08 

842.04 

842.93 

841.40 

1.75 

871.04 

872.11 

873.17 

873.99 

872.42 

2.00 

900.10 

901.10 

902.23 

903.00 

901.41 

3.00 

1000.39 

1001.21 

1003.56 

1004.18 

1002.48 

4.00 

1080. 22 

1080.90 

1083.37 

1083.86 

1082.23 

5.00 

1143.93 

1144.50 

1146.93 

1147.32 

1145.81 

6.00 

1194.83 

1195.29 

1197.61 

1197.92 

1196.53 

7.00 

1235. 52 

1235.90 

1238.05 

1238.29 

1237.04 

8.00 

1268.08 

1268.39 

1270.35 

1270.55 

1269.43 

HI 

(A) 

(C) 

(A) 

(B) 

- ... (c) _ 

CPU/I-0 

(min) 

0.125/1.170 

0.189/1.133 

0.154/1.309 

0.177/1.368 

0.198/1.256 
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TABLE 6 


A COMPARISON OF TEMPERATURE RESULTS 
FOR THE CASES WITH h = h^) 



Case No. 

Time 

(sec) 

12 

13 

16 

17 

18 

0 

293.00 

293.00 

293.00 





0.25 

401.14 

401.27 

401.06 





0.50 

445.89 

446.07 

445.82 





0.75 

478.54 

478.76 

478.61 





1.00 

505.58 

505.83 

505.82 





1.25 

529.69 

530.07 

530.17 





1.50 

552.17 

552.62 

552.80 





1.75 

573.65 

574.15 

574.40 





2.00 

594.44 

594.99 

595.28 





3.00 

673.30 

673.96 

674.38 





4.00 

746.00 

746.73 

747.09 . 





5.00 

813.37 

814.15 

814.46 





6.00 

875.67 

876.52 

876.80 





7.00 

929.48 

930.38 

930.66 





8.00 

977.03 

977.96 

.978.18 





At-sec 

(C) 

(B) 

(B) 

(B) 

(B) 

CPU/I-O 

(min) 

0.224/1.487 

0.200/1.606 

0. 192/1.615 

- 


- 

Parameter P 

0.55. 

0.55 

1.0 

0.25 

0 


134 


























TABLE 7 


A COMPARISON OF TEMPERATURE RESULTS 
FOR THE CASES WITH h = h(T ave ) 


ec; 


Case No. 


14 


15 


0 

293.00 

293.00 

0.25 

350.94 

351.29 

0.50 

387.00 

387.63 

0. 75 

420.17 

420.99 

1.00 

453.04 

454.05 

1.25 

486.41 

488.16 

1.50 

521.83 

524.19 

1.75 

559.79 

562.77 

2.00 

600.52 

603.20 

3.00 

724.14 

725.32 

4.00 

832.04 

833.01 

5.00 

928.13 

928.48 

6.00 

1006.74 

1006.81 

7.00 

1075.18 

1075.07 


1135.05 

1134.79 

At-sec 

(B) 

(C) i 

CPU/I-0 

(min) 

0.192/1.661 

0.235/1.490 
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Figure 1 -vThe Composite Slab '(Pi-P.e.) and its. Finite . Element Model, 

■ * v.-.-vr - . nit ; w-yt i‘\ ■ vx}.i it- ■ k »r,HS- 






$ EXECUTIVE CONTROL DECK 

ID H P LEE, NASA-GSFC, SPACE TECHNOLOGY DIVISION 
A HP Ht«T 

Sul 59 

CEND 

$ CASE CUNTROL DECK 
EC HO= OUT H 
L I NtS=75 

T I TL E= DEMONSTRATION OF TEMPERATURE-DEPENDENT CONVECTIVE F IL M COE FF 1 C I ENT 

SU3 T I Tl E= WITH A CONSTANT HUT GaS TEMPERATURE AT GRID 200 

MPC=303 

DLUAl) = AO 1 

NUiVL INfcAR=A05 

TS T E P = 60 1 

I C= 70 1 

T Emp ( MATE K I AL ) = 705 
UDTPOT 

SET 20=1,2,200 
0L0A0=20 
ELrORC E=20 
NLL0 aD=20 


i>£T *tO — 

1 , A , 7 

,8,100,150 

,200. 


THERMAL: 

= AO 




8EGIN bULK 




S BULK. 

DATA 

DECK 



$ GRID 

POINTS 



GRID 

1 


0.00 


GRID 

2 


0.05 


GRID 

3 


0.10 


GRID 

A 


0.15 


GRID 

5 


0.20 


GRID 

6 


0.25 


GRID 

7 


0.30 


GRID 

8 


0.70 


GRID 

200 


-1.0 


S CONDUCTION 

ELEMENT S 



CRUO 

1 

101 

1 

2 2 

CROO 

3 

101 

3 

A A 

CROD 

5 

101 

5 

6 6 

PROD 

101 

102 

1.0 


MATA 

102 

0.16258 

A. 029 


CROD 

7 

103 

7 

8 

PROD 

103 

10A 

1.0 


MATA 

10 A 

0.00502 

1.725 


S INITIAL CONDITIONS 

1A22. 0 


TEMP 

70 1 

200 


TEMPD 

701 

293.0 


$ ESTIMATED 

TEMPERATURES 


TEMPO 

70 5 

1A22.0 



$ SPECIFIED 

CONSTANT TEMPERATURE FOR THE HOT 

CELAS2 

301 

1.E6 

200 

1 

DARE A 

501 

200 

1 

1A22.E6 

TLUAD2 

AO 1 

501 

0 

o 

• 

o 

+ LD1 

0.0 

0.0 



S TIME 

INCREMENTS AND 

INTEGRATION STEPS 

TSTEP 

601 

25 

0.005 

10 

+ TS1 


35 

0.025 

5 

+ TS2 


20 

0.05 

5 

+ TS3 


120 

0. 10 

10 


101 

101 

101 


3 

5 

7 


1.E5 


0-0 . 


0.0 


+ LD1 

+ TS1 
+ TS2 
+T S3 


enddata 

/- 


Figure 2 — A Listing of the Input Data Deck of the Demonstrative Problem 


137 



S CUNVtOiIVE BOUNDARY IN CONTACT WITH HOT GAS 
CHbDY 1001 1002 PUINT 1 

+CY 200 

PHbUY 1002 1003 1.0 

i-1 AT 4 10 0 3 0.467 

(a) A constant h by the conventional modeling technique 


■ N 0 L INI 40 5 1 1 — i.o 1 

TABGED1 4051 

+TB1 293.0 136.831 589.0 275.063 867.0 

+TB2 1422.0 664.074 £NDT 

’NOG INI 405 1 1 664.074 1 

TA8GED1 4052 

+TB52 293.0 1.0 1450.0 1.0 ENDT 


1 4051 

+ 181 

404.889 1144.0 534.248 +T82 

1 40 52 

+ TB52 


(b) A constant h by the equivalent nonlinear thermal load approach 


N0L INI 

40 5 

1 

1 

1*0 

150 

1 

40 54 



TA'BLEOl 

4054 








+ TB7 

+ TB7 

293.0 

187.704 

589.0 

328.482 

867.0 

446 . 508 

1144.0 

558 .846 

+ Tbo 

+ TB8 

1422.0 

664.074 

ENDT 







NGL INI 

40 5 

1 

1 

-1.0 

150 

1 

40 55 



TABLE01 

40 5 5 








+ T B55 

+ TB55 

293.0 

38.676 

589.0 

136.059 

867.0 

272.238 

1144.0 

449 .59 2 

+ 1 Bt>6 

+ TB 86 

1422.0 

664.074 

ENDT 








(c) h 

= h(Tj_) by the equivalent nonlinear 

thermal load approach 




6 USING 

MFC TO 

CREATE . 

AN AVERAGE ' TEMPERATURE AT GRIU 

100 FOR 

GRIDS 1 AND 200 

MFC 

305 

100 

1 

1.0 1 1 

-0.5 

+ MC 

+ MC 


200 

1 

-0.5 



S CREATING GRIO 

i 150 TU 

AVOID 

GRIO 100 BEING D I KECTG Y 

REFERENCED BY N0GIN1 

GRID 

100 


-0.5 

0.0 0.0 



GRID 

150 


-0.5 

1.0 0.0 



CR0.D 

8 

105 

100 

15o 



PRUD 

105 

106 

1.0 




MAT 4 

106 

5.0 

0.1 




T EH F 

79 1 

100 

857.5 

150 857.5 



c V. . O. st. V. 

4> T "V V 'I s "V V 

v v Vr'M'i'r 

s;c ^ s£ jjs :;s : 


v V v *j! v v W * 5 * W ¥ W "T 'c *r T ^ v 'r v v v ^ 


& SIMUG 

aTING te 

MFERA TURE-UEPENU ENT FIGM C OEEE i C I ENTS 

bASEQ ON 

T ( 150) 

NOG INI 

40 5 

1 

1 

. 1 .00 lbo 1 

406 1 


TABLEU1 

4061 





+ TB1 

+ TB 1 

857 .0 

187.704 

1005.5 

328.482 1144.5. 446.308 1283.0 

. 558.846 +TB2 

+ TB2 

1422 .0 

664.074 

ENDT 


. 


NOG INI 

405 

1 

1 

‘-l.o’ iso . 1 

4062 


TABGED1 

4062 





+ TB3 

+ TtJ 3 

857.0 

113.19 

1005.5 

136. 0S9 1144.5 272.238 1283.0 

449 .59 2 +T 84 

+ TB4 

1422 .0 

667 • 0 74 

ENDT 





(d) h = h(T r ) by the equivalent nonlinear thermal load approach 


Figure 3 - The Images of Packets of Data Cards for Simulating 
Various Convective Film Coefficients 
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$ . T IME 

INCREMENTS 

T S T'frPa 

601 

40 

+*TS3- 


20 

+TS4 


60 


ND INTEGRATION STEPS 
0.025 10 

0.05 5 

0.10 10 


At-set (a) 


+ T S3 
+ TS4 


$ TIME INCREMENTS ANO INTEGRATION STEPS 


TSTEP 

601 

25 

0.005 

10 

; + T S 1 


35 

0.025 

5 

+ T S'2 “ 


20 

0.05 

5 

+ TS3 


120 

0.10 

10 


a t-set (b) 


+TS1 
+ TS2 
+T S3 


$ TIME INCREMENTS IN 
JS.TEP 601 400 

•+T1 120 


INTEGRATIONS AND PRINT-OUT 
0.005 5u 
0.05 20 


a t-set (c) 


+ T1 


Figure 4 - Three Sets of Time-Step Size Used in Computer Executions 



